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The characterization of hot quasi-projectiles produced in symmetric 
or quasi-symmetric reactions (Au + Au, Xe + Sn, Ni + Ni, Ar + KC1) 
at different incident energies are estimated by means of two different 
procedures. The advantages and disadvantages of each method are 
analyzed on the basis of simulations using events produced by two 
slightly different models: HIPSE and ELIE. 



1 Introduction 

Nuclear reactions around the Fermi energy produce hot species by means 
of energy and/or mass transfers between the two partners of the collision. 
Of particular interest is the production of hot projectile-like nuclei because, 
due to favourable kinematical conditions, they can be accurately detected 
by large multidetectors such as INDRA pQ. The thermodynamics of such 
hot species is a key topic of nuclear physics since it allows to address, for 
instance, the existence of a liquid-gas phase transition in nuclear matter. 
The precise measurement of the amount of excitation energy deposited in 
hot nuclei (calorimetry) is a prerequisite in such studies. The complexity 
of nuclear reactions makes such measurements rather difficult and requires 



accurate and well-controlled methods. These latter use event selections and 
kinematical cuts based on the angular distributions of all detected charged 
particles. The aim of the present work is to discuss optimal procedures 
for nuclear calorimetry. For this, computer simulations modeling nuclear 
reactions can be helpfull. Here, we have used the events generated by two 
models: HIPSE (Heavy Ion Phase Space Exploration) [2] and ELIE [3]. 

2 Brief description of the events generators and 
comparison with selected data 

The two event generators are based on a two-step scenario of the reaction: 

• an entrance channel phase ending with the formation (for finite impact 
parameters) of excited PLF (projectile-like fragment), TLF (target-like 
fragment) and participants. 

• a second phase considering secondary decay and propagation towards 
the detectors. 

In the ELIE model, the mass numbers of the PLF, TLF and partici- 
pants are obtained by considering the geometrical overlap of the nuclei for 
each impact parameter according to the high energy participant-spectator 
picture while in the HIPSE case, it is governed by the minimum distance of 
approach obtained with a trajectory calculation based on a realistic inter- 
action potential. 

To build the kinematics of the projectile-like, the target-like and the 
partition of the participants, the following hypothesis are assumed: 

1. The momentum distribution of the incoming nucleons inside the two 
partners is supposed to have no time to relax on a time scale com- 
parable with the reaction time. This is a frozen approximation: only 
a few hard nucleon-nucleon collisions can occur and those latter are 
governed by a single parameter: the mean free path. 

2. In the ELIE model, the partition of the participants is generated by 
a random process in momentum space. The mass number A of each 
species (including A=l free nucleons) is sequentially chosen at random 
by picking A nucleons from the nucleon momentum distribution. For 
IMF's (A > 4), the excitation energy, E*, is obtained by summing the 
center-of-mass kinetic energy of all nucleons belonging to the fragment. 
If E* is larger than a maximum value associated with a maximum 
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Figure 1: Particle multiplicities cross-sections (from protons to IMF as indicated in 
each panel) for Xe + Sn central collisions at 50 A.MeV. Black points are INDRA 
data while the histograms are for ELIE data (Colour online). 



temperature T max , the fragment is rejected and a new try is made 
until all nucleons have been assigned. In the following, it turns out 
that a value of T max = 5.5 MeV allows to reproduce the experimental 
data (with a level density parameter equal to A/10, E*= 3 MeV/u). 

3. In the HIPSE approach, an aggregation algorithm based on cuts in 
momentum and real space is applied to build the fragments and con- 
stitute the partition. 

In a second step, the partition is propagated in space-time and secondary 
decays are considered. This is done using a standard evaporation code: the 
SIMON event generator [3]. The essential difference between the two models 
lies in the excitation energy distribution in the fragments since in the ELIE 
case, the limiting temperature (5.5 MeV) limits E* to around 3 MeV/u while 
for HIPSE, the amount of excitation energy can extend to much larger values 
(see blue histograms in figj3]) (colour online). 



Figure [T] shows a comparison of the multiplicity cross-sections from pro- 
tons to IMF's between INDRA data and ELIE data for Xe+Sn central col- 
lisions at 50 A.MeV. For that purpose, a selection of the events has been 
made using the total light charged particle transverse energy, Ef r i2, nor- 
malized to the CM available energy. Assuming that this latter is correlated 
to the impact parameter, central collisions have been defined by keeping 
only the 50 mbarn cross-section associated with the largest Et r i2 values. 
For such central collisions, particle multiplicities are correctly reproduced 
by the model. Other observables have been considered and a general good 
agreement has been reached between INDRA data and ELIE data. Same 
conclusions can be reached as far as HIPSE is concerned. These models are 
therefore sufficiently realistic to allow a check of the calorimetric methods 
described in the following. 

3 Data analysis and calorimetries 

PLF's characteristics (mass, charge and excitation energy) are measured by 
considering fragments and light charged particles on an event by event basis. 
The question of the undetected neutrons is discussed later. A completion 
criterium is first applied to the data by imposing that most of the charged 
particles (around 80 percent at the forward of the center of mass) emitted in 
the reaction be detected. This ensures that no major biases are introduced 
in the final retained event sample. Then, two different methods are used: 
a standard calorimetry and the so-called 3D calorimetry. The projectile- 
like frame [5j in which all kinematical quantities are calculated is built by 
adding all fragments (charge larger than 2) which are located in the forward 
hemisphere of the center of mass of the reaction (for symetric systems this 
corresponds to positive CM velocities). Considering the standard calorime- 
try, all particles emitted in the forward hemisphere of the PLF frame are 
assumed to be emitted by the PLF and as such are taken into account in the 
reconstruction of the excitation energy and mass. In order to cover the whole 
47r-space, this contribution is doubled on an event by event basis. Note that 
in this case, no attempt is made to consider non-equilibrium effects. 

Such effects are tentitatively taken into account in the 3D calorimetry 
by first assuming that PLF evaporation is cylindrically symmetric around 
the perpendicular axis of the reaction plane, as shown in figEJ It is then 
assumed that all particles emitted in the angular domain (p [0°, —60°] result 
from the evaporation of the PLF, and, as such, the evaporation probability 
in this angular range prob re f(i) of species i (proton, deuteron, triton, ...) 
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Figure 2: Angular domains with respect to the reaction plane and beam direction 
chosen to determine the probability of emission of particles by the PLF. 



is equal to 1 whatever the kinetic energy (E ci ) distribution, called W re f(i). 
In other angular regions covering the same solid angle, other processes of 
emission can occur so that the evaporation probability can take values lower 
than 1. The probability Pj(i) in this region labelled j is then obtained by 
dividing the kinetic energy distribution Wj(i) by W re f(i) for each energy bin 
E c (i). Thus, with this method, the probability of evaporation is a function 
of the species (proton, deuteron, triton, ...), the kinetic energy and the angle 
of emission of the particle. The probability of evaporation hereafter labelled 
probi for particle labelled i is used for the reconstruction of the hot PLF : 
in mass (Aplf), charge (Zplf), and excitation energy E* plf according to 
the following equations: 

multot 



Aplf = Yl P robi X Ai + N neutron (1) 

i=l 

multot 

ZpLF = Yl P r ° bi X Zi ( 2 ) 



i=l 
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Figure 3: < E* PLF > (MeV/u) as a function oi E trl2 for Ni+Ni at 52 A.MeV. Left, 
ELIE data. Right: HIPSE data, both compared with INDRA data (red crosses and 
open points) (Colour online). 
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The number of neutrons is obtained by imposing the conservation of 
the (A/Z) of the projectile. The mass number of fragments when not 
measured is obtained from the measured atomic number following the 
Charity prescription [6]. 



Figure [3] shows preliminary results concerning the correlation between 
Et r i2 and the mean PLF excitation energy per nucleon — for ELIE data 
(left) and HIPSE data (right). In both cases, the INDRA data have been 
superimposed for a sake of comparison. 

Results of the studies are rather contrasted. None of the methods can 
reproduce the excitation energy values given by the ELIE model and a strong 
over-estimation is obtained especially for the most central collisions. This is 
due to the fact that, in the ELIE case, even in the reference angular domain, 
most of emitted particles do not originate from a slow equilibrated process. 
However, for excitation energies lower than 2 MeV/u, the 3D method gives 
good results. As for HIPSE, the 3D calorimetry reproduces the excitation 
energy distribution while the standard method leads to a systematic over- 
estimation. It is worth noted that, whatever the method, the methods 
applied to INDRA data lead to very high excitation energies: a point which 



remains to be understood. 

4 Summary and perspectives 

Two different methods for the characterization of projectile-like fragments 
in nuclear reactions have been discussed. Tests of the procedures using data 
generated by two models, ELIE and HIPSE, have been shown. Results are 
contrasted. However, work is still in progress to improve the 3D calorimetry. 
This could help extract, by subtraction of the evaporation component, the 
non-equilibrium processes associated with the first instants of the reaction. 
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